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Classical dynamics is formulated as a Hamiltoniari flow on phase space, while quantum mechanics 
is formulated as a unitary dynamics in Hilbert space. These different formulations have made it 
difficult to directly compare quantum and classical nonlinear dynamics. Previous solutions have 
focussed on computing quantities associated with a statistical ensemble such as variance or entropy, 
ffowever a more direct comparison would compare classical predictions to the quantum for contin- 
uous simultaneous measurement of position and momentum of a single system. In this paper we 
(•~*) , give a theory of such measurement and show that chaotic behaviour in classical systems can be 

• reproduced by continuously measured quantum systems. 
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I. INTRODUCTION 



The Hamiltonian formulation of classical mechanics assigns simultaneous, arbitrarily accurate, values for the canon- 
ically conjugate position and momentum to distinguishable particles. Indeed in classical mechanics these simultaneous 
values are regarded as properties of the particles themselves; measurement simply reveals these values and need not, 
in principle, add uncertainty to their determination. In quantum mechanics the conjugate position and momentum 
are represented by noncommuting operators, x and p, with [x,p\ ~ ih. It follows that there is no physical state for 
QQ ' which position and momentum can take dispersion free values. This is the content of the standard formulation of the 
, uncertainty principle {(Ax)^) {(Ap)"^) > h^/i, where AA — A— (A). This does not mean, however, that simultaneous 

■ measurements of position and momentum are impossible, only that such simultaneous measurements cannot be made 
', arbitrarily accurate [Q. 

■ Our objective is to show that the observed phase space, reconstructed from the classical stochastic measurement 
C^h' record of continuous joint measurements of position and momentum, corresponds to the classical phase space de- 
^Jj ■ scription with added noise. By a careful specification of the measurement model, including the relevant states of the 
^ ■ apparatus, together with a Markov assumption, we show that the observed classical stochastic measurement record 

corresponds to the conditional quantum averages {x{t)) and {p(t)). We obtain a stochastic Schrodinger equation for 
the conditional state of the system with complex noise, together with classical stochastic equations for the variables 
corresponding to the observed measurement records. This extends the results of Bhattacharya et al. ||^ where only 
a position measurement was considered. Other authors 1^, U, ||, ^, 0, || have considered a general type of stochastic 
Schrodinger equation (which may be interpreted as corresponding to a measurement) in the framework of quantum 
state diffusion to examine dissipative chaos in open systems. 

A model for the simultaneous measurement of position and momentum was given long ago by Arthurs and Kelly 
1^. This model consists of two meters which are allowed to interact instantaneously with the system. The interaction 
couples one of the meters to position and the other to momentum, encoding the results of the measurement in the final 
states of the meters. Projective measurements are then made on each of the meter states separately. In more general 
terms the Arthurs-Kelly model is an explicit example of the Gelfand-Naimark-Segal theorem which enables us to 
write a positive joint probability density of two real variables (the actual measurement outcomes) in terms of an inner 
product on an extended Hilbert space (extended to include the two apparatus as well as the system to be measured) . It 
can also be considered as an example of a positive operator valued measure ^ . The Arthurs-Kelly model suggests 
a route to a phase space description of quantum mechanics based on the joint outcomes of simultaneous measurements 
of position and momentum. However, to achieve a phase space description of the dynamics it is necessary to consider 
a time sequence of such measurements. The Arthurs-Kelly model forces the conditional state of the system into a 
coherent state after a single measurement and is too strong for this purpose. 
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We generalize the Arthurs-Kelly model to allow for a weakening of the measurement. In this way repeated mea- 
surements can take place without invariably reducing the system to a coherent state. By simultaneously increasing 
the number of measurements and weakening the strength of each measurement, the continuous limit is achieved. We 
show that under continuous observation the evolution of the system state is described by an Ito stochastic Schrodinger 
equation. Unlike evolution via the ordinary Schrodinger equation, continuous observation forces the conditional state 
of the system to remain localized. Hence the quantum mean of the phase space variables can be thought of as a 
trajectory. It is in this light that we investigate quantum nonlinear dynamics under continuous observation. 

The theory of continuous measurement in quantum mechanics was first developed by Barchielli and co-workers 
0, [| (see also @ 0) and some of the results in this paper are contained in their work. However our 
approach uses the Arthurs-Kelly measurement model as a starting point and follows the work on continuous position 
measurements by Caves and one of us in In this way our results have a clear physical interpretation. 

In section II we generalize the Arthurs-Kelly measurement model. In section III we derive the stochastic Schrodinger 
equation for continuous measurement which is then investigated analytically in section IV and simulated numerically 
in section V for a chaotic system. Finally, in section VI we discuss our results. 



II. PHASE SPACE MEASUREMENTS 



Following Braunstein et al. ]20| we rederive and generalize the result implicit in the work of Arthurs and Kelly 
. The simultaneous measurement of position and momentum of a one dimensional quantum system is achieved by 
constructing an interaction governed by the Hamiltonian 

HT^H{x,p) + {^^xpi + spp2)d{t~tr) (1) 

which couples the system Hamiltonian H, with two detectors di and c?2. The measurement model also requires that 
the detectors have been prepared in the initial states 

{x,\d.) = (^nA0-^exp(^2^) (2) 
ipM) - f?V%xpf^) (3) 



where Ai = ^ and A2 = The parameter s is called the squeezing parameter and biases the coupling so that 
one may obtain more information on either position or momentum. The parameter a will be used to weaken the 
measurement, decreasing the amount of information collected on both position and momentum. 
Before the interaction the system state p and the combined detector state 

Pd = \dld2){did2\ = ® \d2){d2\ (4) 

are assumed to be uncorrelated with an initial combined density operator oi p® pd- At t = the evolution operator 
for the interaction 

Ui = cxp [--i [\xpi + SPP2)] (5) 

couples the system to the measurement apparatus, entangling the system and detector states. After the interaction 
the probability of finding the detector positions in the small area [xi, xi -\- dxi] x [x2, X2 + dx2] is 



Prob (xi, X2)(ixi(ia;2 — tTtTd[UippdU}\xiX2){xiX2\jdxidx2 (6) 

= tv (t^{xi,X2)p%{xi,X2)''^ dxidx2 (7) 

where tr(. . .) and tTd{. ■ ■) are the traces over the system and detector states, respectively, Prob(xi,a;2) is the 
probability density, and the resolution operator 

'T„ixi,X2) = {xiX2\Ui\did2) (8) 

= {2Trfi)~^ J dpidp2exp[-j-(pi(^x- xi) +p2isp- X2))] {piP2\did2) (9) 

= Disxi, ^X2)f,{0,0)D{sxi^X2V . (10) 



The displacement operator 

D{fi, v) = cxp [-j- {jip - vx)] 

and 

t40,0) = {27rh)-^ J dpidp2{pi,-p2\did2)I){sp2, ^^PiV . 
We now define the annihilation operator 

^=^(^^^ + ^75^) 

which satisfies 

[a, at] = 1 

and rewrite the displacement operator as 

D{C) = exp (Cat „ 
where C — (^^/sfi + • It can be easily shown that 

^(C)t = ^(-C) = ^(C)-^ 

D{C)aD{Cy =a-C. 

The coherent states are now defined by applying the displacement operator onto the vacuum 

\a)^D{a)\0). 

They satisfy the following relations 

a\a) = a\a) 
Di(3)\a) ^exp[^{a*f3-a(3*)]\a + p) 
(a|/3) ==exp(-i|a|2-i|/?|2 + a*/3) 



TT-^ j d^a = 

Now defining g = {^\/sp2 + can rewrite ( |l^ ) as 

t„(0,0) = {2T:h)-ia7r-^ J d^ g{ag\0)D{g)'' . 
Using the above relations and integrating we obtain 

(a|t. (0,0)1/3) = (27r?i)-^^exp(-i|a|2-i|/3|2 + £!^a*/3' 

= (a|(27r?i)-^^ :exp(-^at, 
where : • : denotes normal ordering. Hence, 

t<,(0,0) = {27rh)-i^ : exp (-^dh 



If we now define X ^ Xi + *X2 = (^^/sxi + i-^X2^ then the resolution operator becomes 



= (2nr"^TAx)- 
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The probability of finding the detector positions in the small area [xi, Xi + d-Xi] x [X2i X2 + dx2\ is now 

PYoh{x)d^X = 2nPmh{xi,X2)(fx (31) 

= tr (t,{x)p%{x)^) d\ (32) 

= tr(4(x)/5)d'x (33) 



where rf^x = c^X 1*^X2 and 



4(X) = t,(x)^T,(x) (34) 



-^(^)':exp (^)^(at-x*)(a-x) 

is an effect density |p4| . It can be easily shown that 

x"F,(x)d2x = (36) 

lxP4(x)rf\ = ««^ + (47^)'- (37) 
Hence, defining the notion of a mean for this measurement process 

{f{x))a = I /(x)Prob(x)d\ (38) 

/(X)tr(4(x)p)rf'x (39) 
we find that 

(X)^ = {a) (40) 

(IxP). = (aat) + (47^)' (41) 

or 

(2:1). = (X), (X2)a = (p) (42) 

= (43) 

(X2^).-(P^> (44) 

where (A) = tr(A/5) is the quantum expectation. Thus the readout variables xi and X2 give, respectively, the position 

and momentum of the system with additional noise dependent on cr. We obtain maximal information from the system 
when the variances are at a minimum. That is, when a = 1. In this case 

Fi(x) = TT-i : exp [-(at - x*)(a- x)] : (45) 

= vr-i|x)(xl (46) 
and the probability density reduces to the Husimi density or Q function 

Prob(x) = 7r-i(x|/5|x) (47) 

= 7r-iQ(x). (48) 

Suppose we take a measurement and obtain the outcome x'- As a consequence of the strength of this measurement 
the system state collapses to 

P' - l^^^(x')ptr(x')^ (49) 

= r^^lx')(x'l/5|x')(x'l (50) 
(x IpIx ) 

= lx')(x'|. (51) 



(35) 
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However when cr ^ 1 the resohition operator has the asymptotic expansion 



T.(x) 



TT 2 . 



-^(l + 2(at-x*)(a-x))J + O(^). (52) 
Using this result one can show that the system state conditioned on the measurement outcome x' has the expansion 
p' = p+ ^{ tr [(at - - X')p] - («t - x'*)(« - X'), P} + O(^) (53) 

and thus, if a is large enough, the process of measurement will have negligible effect on the system. 



III. CONTINUOUS MEASUREMENT 



Consider a sequence of phase space measurements governed by the Hamiltonian 



H{x,p) + {-^xpi + SPP2) ^ '5(t - nSt) 

n=0 



(54) 



where after each measurement the detectors are reset into the initial states given by equation or (^). The 
assumption that the detectors are reset is equivalent to making a Markov assumption for a single apparatus coupled 
to the system. By resetting the detector at each time step we ensure that no coherent memory of the system state 
survives in the states of the apparatus. Following we will first derive the master equation for unconditional (or 
nonselective) evolution of the system density operator in the continuous limit i5t 0, cr —> 00. By unconditional 
evolution we mean that no account is taken of the measured results. Thus after each measurement occurs we ignore 
the result and average over all possible measurement outcomes. If we denote the system density operator immediately 
before the n-th measurement by p{n5t) then 



p{n5t + 5t) = U f d^x%{x)p{nSt)%{x)^m 



where 



= 1 - j-H6t + 0{5t^). 



For any operator A it is possible to show that 



A 



(55) 

(56) 
(57) 

(58) 



Hence we obtain 

p{n6t + 5t) - p{nSt) 
St 



By setting t — nSt and taking the continuous limit 5t ^ 0, a 
master equation for unconditional evolution: 



H,pinSt)\ ~ [a, {a\p(n8t)\\ + 0{8t) + 0( J^) + 0(3^). (59) 
00, with 7 = l/(5t(7^ held constant, we obtain the 



dp 
'dt 



H,p 
H,p 



7 [a, 



x,p 



2h 



p, 



(60) 
(61) 



where Fi = 7s and r2 = ^. This equation has already been derived by Barchielli et al. ||l3|, but their approach is 
somewhat different. By setting F2 — in (61) we obtain the unconditional master equation for continuous position 
measurements previously derived in [p^ . 

We now wish to derive the conditional (or selective) master equation for the system density operator. In this case 
the evolution of the system is conditioned on a history of measurement readouts 



{x{n6t)} = {x{0),x{St),x{m,---} 



(62) 
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where x(n5t) is the detector position for the n-th measurement. Hence, if p{nSt) is the system density operator 
immediately before the n-th measurement, then 

p{nSt + 6t) = U[Pmh{x{n6t))]'^%{x{nSt))p{nSt)%{x{nSt))^Ul (63) 
To proceed we extend the definition of the readout variable by setting 

X{t) = xinSt) for nSt<t <{n + l)6t (64) 

and introduce the new variable 

Xit)= f x{t')dt'. (65) 
Jo 

Hence A'(O) = and 

n-l 

X{n5t) ^5t^ x{mSt) for n > 1. (66) 



Using (|^ and (El|) we obtain 



E,{x{n5t)) = {x{n5t))<r (67) 

= tr{ap{n5t)) (68) 

V,{x{n5t)) = {\x{n5t)~{x{n5t)),\^)^ (69) 

= {\x{n5tr).-\{x{n5t))^^ (70) 

{aa^p{n8t)) \ tr (ap(nJt))|' + (4^) ' . (71) 



tr I 



where the subscript c has been added to emphasize that the mean and variance are conditioned through p on the 
entire history of measurement readouts. Now letting 

5X{n5t) = X{n5t + 5t) - X{n5t) = 5tx{n5t) (72) 

we find that 

Ec{6X{n5t)) = 5tiv{ap{n6t)) (73) 

Vc{5X{n5t)) = (^tr(aatp(„5t)) _ |tr(ap(ri5<))|V (4?^)^^ (74) 

= \-i-^5t + 0{5e). (75) 

where we have set 7 = l/Sta"^ to be constant in anticipation of the continuous limit. Hence for 5t small enough we 
can approximate increments in the variable X by 

5X{n5t) = tr {ap{n5t)) St + ^j-i5^{n6t) (76) 

where it is understood that the complex Ito increment d£, is of order St^ and satisfies E (5^) ~ 0, E {S^*S^) = St. In 
the continuous limit St —> with t = nSt constant we have 

dX{t) = tr {ap{t)) dt + ^-f-^d^{t) (77) 

where ^(t) is a complex Wiener process and the Ito differential satisfies the algebra 

E{d() = (78) 

dCd^ = dt (79) 

= (80) 

d^dt ^ 0. (81) 
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To simplify the following we will always replace S£^*S£, by 5t and set = d£,*^ = in anticipation of the above 
algebra in the continuous limit. Using (72) together with ( |76|) we obtain the following asymptotic expansion for the 
resolution operator (pO) 



T.(x) 



( 2^ 



\ 7re 



(l + 7^ (i^e + ASC) - l5t (iU + \))+ 0{St^ 
where 

A = a — tr (ap) 

and it is understood that p — p{nSt), x — ^^J^d ~ S^{n6t). Hence we find that 



%ix)p%{xf = 4^ (p + 7* {aU^ + ASC,p} ~ l5t 
and thus, using (|57|) we obtain 



it P 



-4, 



it 



H,p 



+ {aU^ + ASC, p} - jSt 
~j6t [a, [at,p]] +jiH[a'']pd^ + jiH[a]pSC + 0{dti) 



where we have defined the superoperator 

n[A]p={A^ tr(ip),p}. 

In the limit St ^ we obtain the master equation for conditional evolution: 
dp{t) = p{t + dt) - p{t) 

H, pit)] dt~-i [a, [d\p{t)]\ dt + -i^rL[d^]p{t)di{t) + -i^'H[d]p{t)di{tY 



(82 
(83 
(84 

(85 
(86 

(87 



(89 



It is easy to see that upon averaging this stochastic differential equation we reproduce our original master equation 
for unconditional evolution (|60|). However, note that unlike the unconditional equation, this equation preserves the 
pure state property of p. One can easily prove this by showing dp^ = dp under the assumption that p = 
where dpp' = {dp, p} + dpdp. As a consequence of this fact, the above master equation has an analogue for pure state 
evolution in terms of a stochastic Schrodinger equation: 



rfl^W) = -TiH\i'{t))dt - 7 (fltfl + 1 _ (at)a _ a) {d) + \{d)\^) \i^[t))dt 

+7^ (at - (at)) \ij{t))dm + 7^ (a - (a>) mmit)' 



(90) 



where {A) — {ij}{t)\A\^{t)) . In terms of position and momentum variables the master equation for conditional evolution 
reads as 



dp 



H,p 



dt~^T, 



x,p 



dt~^T, 



P,P 



dt + n^^Ti^n[x]pdWi + h^^V2^n\p\pdW2 



n) 



where 



E{dWi) = 
dWidWj = Sijdt 
■^{dWx+idW2) = d^ 

and the readout variables xi and X2 obey the following stochastic processes 

dXi{t)^ tv{xp{t))dt+ \h^Ti-^dWi{t) 
dX2{t) = tr {pp{t)) dt+ ^hh2-^dW2{t) 

with 



Xi{t)= f Xi{t')dt' X2{t)= f X2{t')dt'. 

Jo Jo 



(92) 
(93) 
(94) 



(95) 
(96) 

(97) 
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By setting r2 = in ( pl\) we obtain the conditional master equation for continuous position measurements previously 
derived in [l9[| . Note that Xi and X2 are charged by stationary white noise 

xi = {x)+ ia^rrwi (98) 

X2 = {p)+ ^hh2-iW2 (99) 
making their graph highly irregular. It is thus better to represent the measured trajectory by (x) and (p). 



IV. ANALYTICAL INVESTIGATIONS 

We will now investigate the effect of measurement on the system state by setting H — Q and defining 

K = E{{x^)~{x)^) (100) 
V, ^ E{{f)-{p)^) (101) 

C^p = E(^\{xp) + \{px)~{x){p)) (102) 

where Vx and Vp are the expected variances in position and momentum, and Cxp is the expected covariance between 
position and momentum. The average is taken over all possible measurement histories. One can then derive the 
following set of coupled differential equations: 

^ ^ 2^ _ ^Vx' - ^Cxp' (103) 
at s n sn 

— ^ s^n- -Vp - —Cxp (104) 

"^^^^ - -^Cxp{sVx + -Vy. (105) 



dt n '^x s 



The solutions of which are 



Mt) = TTT 7 Z (106) 



^ (2sV° + h tanh(27t)) (^sh - 


f 2V^ tanh(27i)) 


- 4:sC^p^ tanh(27t) 


2« (h + 2sV^ tanh(27i)) {^sh 4 


- 2V^ tanh(27<)) 


-4sC2/ tanh2(27i) 


h 

— as i — > 00 
2s 






(2V° + s?itanh(27t)) (h + 


2sV° tanh(27i)) 


- 4sC°p^ tanh(27t) 


(h + 2sV^ tanh(27t)) (^sh 4 


- 2V^ tanh(27t)) 


-4:sC^p^ tanh2(27i) 


sh 

— ^ — as t —^ 00 
2 






sh C^j 


,sech2(27t) 





^p(t) = VT V Z (107) 



C.pit) - 7 ^7 ^-^ '—^ (108) 

( h + 2sV^ tanh(27i) ](sh + 2V^ tanh(27<) j - 4sC0p^ tanh^ {2-ft) 



—^0 as t 00 

where Vx{0) = V^, Vp{0) = Vp and Cxp{0) — C^p- Hence the process of measurement induces the system state to 
collapse into a coherent state. If the measurement retrieves no information on momentum, i.e. r2 = 0, then 

, > n 4C!?„ Fit 

Vpi^)-K^^r,t-^^^^ (110) 

f,r<0 

C(t) = ^. — (111) 
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and the growth in the momentum variance is unbounded. Similarly, if the measurement retrieves no information on 
position then the position variance grows unbounded. However, when both position and momentum are measured 
simultaneously the system state is forced into a coherent state. When H ^ we expect that for a suitable choice of 
7, any spreading of the quantum wavepacket caused by nonlinearities in the Hamiltonian will be counteracted by the 
measurement induced localization. 

One might naively assume that if the measurement only retrieves information on position (or momentum) then 
the state will not localize. However this is not always the case. Note that when \C^p\ > h/2 in equation ( |110| ) 
the momentum variance will initially decrease. Thus if the system dynamics is such that it increases the covariance 
between position and momentum, then the continuous measurement of position may also localize momentum. For 
example, consider the Hamiltonian describing free particle motion, H — afp' . When r2 = the variances and 
covariance satisfy 

-^K'+4aC,p (112) 

ri?i-^c,p' (113) 
4ri 

^C^pV^ + 2aVp. (114) 

Although we could not solve these equations analytically, it is easy to see that all physical solutions asymptotically 
attract to the stable fixed point 




(115) 



Hence, measurement of position does not introduce a diffusion in momentum, and the state localizes (this result has 
been derived previously in p^). However, for free particle motion, if we only measure momentum the state does not 
localize. The system dynamics accelerates the growth in the position variance. See [§, |5[ 0, ||, |2j, ^ for more 

on localization. 



dt 

dVp 
dt 



dt 



V. NUMERICAL SIMULATIONS 

We will now numerically investigate the solution of the stochastic Schrodinger equation (|9^) for the driven system 

H ^ ap^ + bx^ + ex* + dx cos(wi). (116) 

The numerical method to solve this equation is simple. To take advantage of the measurement induced localization 
we use a local moving number basis 

|n)^^at"|o) n = 0,l,...iV (s = 1). (117) 

truncated at some finite value A^. The stochastic terms are integrated using the first-order Euler method while other 
terms are integrated by diagonalizing the position and momentum operators and using the split-operator formula. 

We will first consider an integrable case when a = b = 5, c=l, d = and Ti = 0.05. The initial state was 
chosen to be a coherent state (s — 1) centered at {x,p) = (—2,1) when t = Q. The Husimi density of the initial 
state together with the contours of the Hamiltonian are plotted in Fig. 1(a). The Husimi density of the evolved 
state {t — 4) together with the trajectories {{x), {p}) for different measurement schemes is plotted in figures l(b)-(d). 
The evolved state in Fig. 1(b) is the result when no measurements occur (7 = 0). In this case, nonlinearities in the 
Hamiltonian cause the state to shear as it evolves, spreading it along the contours. The trajectory has little meaning 
when 7 = 0. In Fig. 1(c) the evolved state is the result of continuous simultaneous measurement of position and 
momentum with 7 = l/\/2 and s = 1. In this case the state has remained localized as it follows the contours. The 
continuous measurement of position only (r2 = 0) when Fi = 1 has also kept the state localized. This is shown in Fig. 
1(d). The combined variance of position and momentum is plotted in Fig. 1(e). We must emphasize that only when 
both position and momentum are measured together does the trajectory correspond via ( p^ , |9^ ) to the outcome of an 
actual measurement. If only position is measured, only (x) observed while (p) is simply the result of a mathematical 
calculation. 

Now consider the chaotic case when a = 5, 6 = —8, c — 1, d — 15, uj = 27t and h = 0.05. A Poincare stroboscopic 
map with unit strobing frequency is plotted in Fig. 2(a). For an initial state the same as above, the evolved state 
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{t = 5) together with the trajectories for different measurement schemes is plotted in figures 2(b)-(d). When no 
measurement occurs (Fig. 2(b)) the chaotic action of stretching and folding spreads the state across the phase space. 
However when we continuously measure position and momentum (Fig. 2(c)) the state remains localized and the 
trajectory resembles classical motion with noise. This noise will vanish as we approach the classical limit h —^ 0. In 
Fig. 3(a) we have plotted the quantum trajectory for the same parameter values as above except with h = 10^®. In 
this case the total variance remains below 10~^ (Fig. 3(b)) and the noise is not visible. The corresponding classical 
trajectory is plotted in grey and is only visible when it deviates from the quantum at t sa 5. The evolution of this 
system under the continuous measurement of position has already been studied by Bhattacharya et al. |^ . They also 
find that the measurement keeps the system state localized. It is not surprising that this is also the case when only 
momentum is measured (Fig. 2(d)). 



VI. CONCLUSION 

We have derived an Ito stochastic Schrodinger equation (Q) describing the evolution of a quantum system under 
the continuous simultaneous measurement of position and momentum. The outcome of this measurement is a classical 
stochastic record obeying (^. As a consequence of the measurement, the system state is forced to remain localized 
allowing a classical interpretation of the quantum mean of the phase space variables as the trajectory of the system. 
This trajectory corresponds to the actual measured trajectory minus noise ( ^ , p9| ). Furthermore, the localization 
property allows a well-defined classical limit via Ehrenfest's theorem. Indeed, for small ?i, numerical results show that 
the quantum system approximately follows classical trajectories. However a more complete theoretical understanding 
of the classical limit under continuous measurement is needed. 
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VII. FIGURE CAPTIONS 



FIG. 1. (a) Contours oi H = -|- bx^ + and the initial state {fi = 0.05). (b-d) The trajectory and final state at 
t = 4 when (b) 7 = 0, (c) 7 = 1/^2 and s = 1, and, (d) Fi = 1 and F2 0. (e) The combined variance of position 
and momentum. All quantities are dimensionless. 
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FIG. 2. (a) Poincare map for H = 5p'^ - 8x'^ + x'^ + 15a; cos(27ri) and the initial state {h = 0.05). (b-d) The 
trajectory and final state at t = 5 when (b) 7 = 0, (c) 7 = l/\/2 and s = 1, and, (d) Fi = and r2 = 1. (e) The 
combined variance of position and momentum. 

FIG. 3. (a) The quantum (black) and classical (grey) trajectories for the same Hamiltonian as in Fig. 2 except 
with h = 10~^. (b)The combined variance of position and momentum. 



